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MULTI-BASE REPRESENTATIONS OF INTEGERS: ASYMPTOTIC 
ENUMERATION AND CENTRAL LIMIT THEOREMS 

DANIEL KRENN, DIMBINAINA RALAIVAOSAONA, AND STEPHAN WAGNER 


Abstract. In a multi-base representation of an integer (in contrast to, for example, the binary 
or decimal representation) the base (or radix) is replaced by products of powers of single bases. 
The resulting numeral system has desirable properties for fast arithmetic. It is usually redundant, 
which means that each integer can have multiple different digit expansions, so the natural question 
for the number of representations arises. In this paper, we provide a general asymptotic formula 
for the number of such multi-base representations of a positive integer n. Moreover, we prove 
central limit theorems for the sum of digits, the Hamming weight (number of non-zero digits, 
which is a measure of efficiency) and the occurrences of a fixed digits in a random representation. 


1. Introduction and Background 

A numeral s?/.sfe7?f] (also called system of numeration ) is a way to represent numbers. The 
most common examples are, of course, the ordinary decimal and binary systems, which represent 
numbers in base 10 and 2, respectively. Besides those “standard” systems, there is an immense 
number of other numeral systems. 

For fast arithmetic, the right choice of numeral system is an important aspect. The algorithms 
we have in mind here are, for example, exponentiation in a finite group and the scalar multiplication 
on elliptic curves. Both are used in cryptography, and clearly we want to improve on the running 
time of those algorithms (which are often based on a Horner scheme, cf. Knuth [l4]). 

Starting with the binary system, one can improve the performance of the aforementioned 
algorithms by adding more digits than needed. Thus, we make the numeral system redundant , 
which means that each element can have a lot of different representations. For instance, using 


digits 0, 1 and —1 can lead to a speed-up, cf. Morain and Olivos 20 for such a scalar multiplication 
algorithm on elliptic curves. To gain back the uniqueness, additional syntax can complement the 
redundant system. In the example using digits 0, 1 and —1, this can be the non-adjacent form, see 
Reitwiesner’s seminal paper 25 . Generalizations in that direction can be found in j3 j9,19 27 . 


A different way to get redundancy, and thereby a better running time of the algorithms mentioned 
above, is to use double-base and multi-base numeral systems. For example, we can represent a 
number by a finite sum of terms ai S ae 7^ e ll le for some digits a^, which leads to a multi-base system 
with three bases. A formal definition is given in the next section. Note that multiplication by one 
of the bases (in the example: 3, 7 or 11) is extremely simple for such representations, just like 
doubling is easy for binary representations. This is a very desirable property for fast arithmetic. 

Double-base numeral systems are used for cryptographic applications, see for example [ljj5,6 . 
The typical bases are 2 and 3. With these bases (and a digit set containing at least 0 and 1), 
each positive integer has a double-base representation, cf. Berthe and Imbert [2|. When using 
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general bases, less is known on the existence, cf. Krenn, Thuswaldner and Ziegler 1161 for some 
results using small symmetric digit sets. However, choosing the digit set large enough (so that the 
numeral system with only one of the bases can already represent all positive integers), existence 
can always be guaranteed. Thus, when each positive integer has a multi-base representation, a 
natural further question arises— and this is also the main question studied in this article: how 
many representations does each integer have? Our Theorem [I] provides an (asymptotic) answer to 
this question. 

The question is also motivated by the cryptoanalysis of evaluation schemes (e.g. elliptic curve 
scalar multiplication): One can avoid side channel attacks if the corresponding numeral system is 
very redundant, i.e., if each element has many different representations. In addition to the number 
of representations, other parameters, such as the (Hamming) weight or the sum of digits, are of 
importance in this context and therefore studied here as well. The Hamming weight in particular 
is a measure for the efficiency of a digit representation for fast arithmetic. We show here that the 
sum of digits and the Hamming weight (as well as the number of occurrences of any fixed digit) of 
a typical representation of n is of order (logn) m , where m is the number of bases. 

Our paper is structured as follows. The following section provides more precise definitions and 
reviews existing results on the number of representations (which are available in very special cases). 
This is followed (in Section [3]) by the precise statements of our main results. These results also 
include, apart from the asymptotic enumeration of multi-base representations, the analysis of the 
sum of digits, the (Hamming) weight and the number of occurrences of a fixed digit. The remaining 
parts of this article (Sections |d] to |8| are devoted to the proofs of all these results, which are based 
on generating functions and the saddle-point method. Section [9] concludes the paper. 

An extended abstract of this paper was presented at the AofA 2014 conference in Paris, see 


15 


2. Terminology and Existing Results 

In a multi-base representation of n (or multi-base expansion), a positive integer n is expressed 
as a finite sum 

L 

n = ^2a £ B £ , O) 

£=1 

such that the following holds. 

• The an (called digits) are taken from a fixed finite digit set T>. Here, we will be using 
the canonical digit set {0,1,..., d — 1} for some fixed integer d > 2, but in principle our 
methods work for other sets as well. 

• The Be are in increasing order (i.e., B\ < B 2 < ■ ■ ■ < Bjf) and taken from the set 

S = {p“ 1 P2 2 -*”:a»€NU{0}}. 

The numbers p±, p m are called the bases (in our setting, these are fixed coprime 
integers greater than 1). The sequence of all elements of S in increasing order is sometimes 
called a Hardy-Littlewood-Polya-sequence. 

In the following, we will discuss the number of representations of n in a given multi-base system, 
which we denote by P{n) (we suppress the dependence on /q, p- 2 , ..., p m and d). Note that this 
number is finite, since our digit set does not contain negative integers. 

For redundant single-base representations a lot is known. Reznick [26 presents results on 
certain partition functions, which correspond to representations with non-negative digits; see also 


Protasov 23 ; 24 for more recent results on the number of representations P(n). When negative 


digits are used as well (for example in elliptic curve cryptography), there are usually infinitely 
many representations of a number, so counting these does not make sense. In this case, expansions 
with minimum number of non-zero digits are of interest, since they lead to fast evaluation schemes. 
See Grabner and Heuberger 10 for a result counting minimal representations (one minimal 


representation is the non-adjacent form mentioned above, cf. also 11 12 25 ). 

Let us consider double-base systems in particular, and let us take bases 2 and p, where p > 1 
is an odd integer, and digits 0 and 1. We can group terms involving the same powers of p and 
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use the uniqueness of the binary expansion to show that double-base representations with bases 2 
and p are in bijection with partitions into powers of p , i.e., representations of the form 

n = ?i,o + mp + n 2 p 2 + n 3 p 3 H- 


with (arbitrary) non-negative integers ne. More generally, the same is true for double-base 
representations with bases q and p and digit set {0,1,..., q — 1}. It seems that the first non-trivial 
approximation of P(n) in this special case is due to Mahler 117 . By studying Mordell’s functional 
equation, he obtained 

log P(pn) ~ (log n) 2 /(2 log p). 

The much more precise result 


log P(pn) = 


2 log p 

-U + 


log,- 

log 71 

log log P 
log p 


1 

log p 


log log n + 0(1) 


log log p\ 

— - log n 

log p ) 


was derived by Pennington 22 . The error term in the previous asymptotic formula exhibits a 
periodic fluctuation. Note that for bases 2 and p , the function P(n) fulfils the recurrence relation 


Pin) 


P(n — 1) + P(n/p) if p | n, 
P(n — 1) otherwise, 


which has been known for a long time in conjunction with partitions of integers. 

For further reference and more information see A005704 in the On-Line Encyclopedia of Integer 
Sequences 21 and see also [5 18 for the connection to double-base systems. 


3. Main Results 

We present our main results now. The aim of this work is to give an asymptotic formula in a 
more general set-up. Throughout this paper, d > 2 and m > 2 are fixed integers, and pi, p 2 , ■ ■ ■, 
p m are integers such that 1 < p\ < p 2 < ■ ■ ■ < p m and gcd(pt,pj) = 1 for i ^ j. As our first main 
theorem, we prove an asymptotic formula for the number of representations of n of the form pi| . 
It will be convenient to use the abbreviation 

__ log d -py 1 

m\ log pj 

Theorem I. If m > 3, then the number P(n) of distinct multi-base representations of n of the 
form (Q satisfies the asymptotic formula 

log P(n) = C 0 (log n) m + Ci (log n)" 1-1 log log n + C 2 (log?z) m “ 1 + C((logn) m_2 loglogn) 
for n —> oo, where 

Co = k , 

Ci = —m(m — 1 )k, 

m 

C 2 = nm(l + \ log pj - \ log d - log (icm)j . 

1=1 

In the case that there are precisely two bases, we have the following more precise asymptotic 
result. 

Theorem II. If m = 2, then the number P(n) of distinct multi-base representations of n of the 
form ([*]) satisfies the asymptotic formula 

P(n) = K(n) (log n) K °n Kl exp( ftlog 2 
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for n —> oo, where K(n) is a fluctuating function of n, that is, bounded above and below by positive 
constants, and 

Kq = \ +2K(log(2At) - 5 (logPi + logp 2 -logd)), 

ICi = 2n(l - log(2«) + ^ (log pi + log P 2 - log d)) - 1. 

Note that the first two terms of the asymptotic formula in Theorem [I] coincide with those in 
Theorem m 

Moreover, we study the distribution of three natural parameters in random multi-base repre¬ 
sentations, namely the sum of digits, i.e. a\ + a 2 + • • • + Ol in the notation of ([^|, the Hamming 
weight (the number of non-zero coefficients af) and the number of occurrences of a fixed digit b. 
We get the following theorems. 

Theorem III. The sum of digits in a random multi-base representation of n of the form ( f^j ) 
asymptotically follows a Gaussian distribution with mean and variance equal to 

hn = (log n) m + 0 ((log n) m ~ 1 log log n) 

and 

ol = ^ 12 log ^ 1 ^ ( ' 1 ° S m + ° (( l0g ^^ 1 l0g l0g 

respectively. 

Theorem IV. The Hamming weight of a random multi-base representation of n of the form ( [¥] ) 
asymptotically follows a Gaussian distribution with mean and variance equal to 

hn = (log n) m + O ((log nr- 1 log log n) 

and 

o"n = ^^j(logn) m + C>((logn) m_1 loglogn) 

respectively. 

Theorem V. Let b £ {0,1,.. .,d — 1}. The number of occurrences of the digit b in a random 
multi-base representation of n of the form (|*]) asymptotically follows a Gaussian distribution with 
mean and variance equal to 

Mn = d d (log n) m + O ((log n) m ~ 1 log log n) 

and 

crl = ^~ g 1 j (l°g») ,n + O ((log n) m ~ l log log n) 

respectively. 

The proofs of all these theorems are based on a saddle-point analysis of the associated generating 
functions. As it turns out, the tail estimates are most challenging, especially in the case m = 2 
(see Section [5] for details). For the asymptotic analysis of the various harmonic sums that occur, 
we apply the classical Mellin transform technique, see 7 . 

4. The Generating Function 

We start with a generating function for our problem. As mentioned earlier, we define the set 

S={p?'p?---p'Zr:a j eNU{0}}, 

which is exactly the monoid that is freely generated by p\, P 2 , • • •, p m - Note that the representations 
of n correspond exactly to partitions of n into elements of S where each term has multiplicity at 
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most d— 1. The generating function for such partitions, where the first variable 2 marks the size n 
and the second variable u marks the sum of digits, can be written as 

f(z,u )=n ( 1 + uz k +u 2 z 2k +• • •+= n (4 t) 

fa£«S h£S 

Likewise, we have the generating function 

G(z,u) = ]^[ (l + uz h + uz<2h H-b uz^ d ~ 1)h ^j = ( 1 + uz h ———^— j , (4.2) 

hes K hes ' ~ « J 

where the second variable marks the Hamming weight (number of non-zero digits, or equivalently 
number of distinct parts in a partition). For a digit b £ {0,1,..., d — 1}, whose occurrences will be 
marked by u, we use the generating function 

H b (z, u) = n (i + zh + • • • + uz bh + • • • + z^ h ) = n ( 1 + (« - l)z bh ) ■ (4.3) 

h£S heS ' ' 

Obviously, F(z, 1) = G(z, 1) = H b (z, 1). We would like to apply the saddle-point method to 
these generating functions. The trickiest part in this regard are the rather technical tail estimates, 
especially when m = 2, which will be discussed in the next section. We will also need an asymptotic 
expansion in the central region. To this end, we define the three functions 

f(t,u) = logF(e~ t ,u) = ^log(l + ue~ ht + u 2 e~ 2ht + --- + u d - 1 e- {d ~ 1)ht ), 

hes 

g(t , u) = log G(e-\ u) = ^ log(l + ue~ ht + ue~ 2ht + • • • + ue~^- 1)ht ) 
hes 

and 


hb(t,u) = logH b (e *, u ) = £log(l + e ht 4 - Vue bht +-he 1 '> ht y 

hes 


Lemma 1 . Suppose that u lies in a fixed bounded positive interval around 1, e.g. u £ [1/2,2]. 

(1) For certain (real-)analytic functions $ 2 ( 0 ), • • •, fm(u ) with 

m . 

fm(u) = log(l + U~\ -h U d ~ 1 ) TT - - , 

jJi l°g p j 

we have the following asymptotic formula as t —> 0 + (t positive and real), uniformly in u: 

m «) = aog i/tr+ /M^ao g w 1 + • • •+/i(«)(io g m +. 

ml (m — lj! 

Moreover, 

= ~ (m"- 1 )! t (1 ° S 1/ ^ )TO ~ 1 + 0(t-\l o gl/tr- 2 ) 

and 

"> = (J-ltiid 108 1/ ' r_1 + °< r2 (log 1/f T~ 2 ) ■ 

Finally, there exists an rj > 0 such that for complex t with |Imi| < rj, we have 
^ 3 /(*.«) = C’((Re<) _ 3 (logl/(Ret)) m ^ 1 ) 
as Ret —> 0 + , again uniformly in u. 
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(2) Likewise, there exist functions gi(u), 52 («), ..., gmip) such that 
9m(M), u\m 1 9m— l( M ) 


ml ym — 1 )! 

and the same conditions as in (1) hold with 


9 i(u)(\ogl/t) +0(1) 


9m(u) = log(l + (d- l)zz) J] 


log Pj 

(3) Moreover, for each digit h £ {0,1,d — 1}, there exist functions h b ,\{u), h b , 2 {u), . 
hb, m {u ) such that 


hb (t,u) = ^M(i 0 gi/*r + ^"^ (logiAr- 1 

ml [m — lj! 

and the same conditions as in (1) hold with 

m 

h b ,m{u) = log (d - 1 + u) 

i=i 


h btl (u)(logl/t) + 0 ( 1 ) 


log Pj 


Proof. To prove the first part, we apply the classical Mellin transform technique to deal with 
the harmonic sums, see the paper of Flajolet, Gourdon and Dumas [7]. Consider first the Mellin 
transform 

f-OO 

Y(s, u) = / log(l + ue-* + u 2 e~ 2t + • • • + u d ~ 1 e~ < ' d ~ 1 ' ,t )t s ~ 1 dt. 

Jo 

Integration by parts allows us to provide a meromorphic continuation (cf. Hwang [13j). We have 


y(s, zi) = 


+ 2u 2 e~ 2t + • • • + (d - 1 )zz d - 1 e-( d - 1 )* 
1 + zze _t + • • • + 


dt, 


which exhibits the pole at 0 with residue log(l + u + u 2 + ■ ■ ■ + u d 1 ), i.e., we have 


T(s, zz) 


log(l 


d -1 


as s —> 0. By repeating this process one obtains a meromorphic continuation with further poles at 

- 1 ,- 2 ,.... 

Moreover, since the integrand in the definition of Y(s,u) decays exponentially as Ret -A oo, we 
can change the path of integration to the ray consisting of all complex numbers f with Arg t = e > 0, 
where e is chosen small enough so that there is no t with Arg t < e for which the expression inside 
the logarithm vanishes (this is possible since u was assumed to be positive, so there are no real 
values of t for which this happens). Set (j = e le , and perform the change of variables t = j3v to 
obtain 

/*oo 

Y(s, u)=P s / log(l + ue~P v + u 2 e~ 2f}v + ■ • • + u d ~ l e -( d - 1 )^)z; s - 1 dv. 

Jo 

If now s = a + ir with a > 0, then the integral is uniformly bounded in r for fixed a, while the 
factor /3 s = e le<T_£T decays exponentially as r —> oo. The same can be done for a = 0 and for 
negative values of <r (after suitable integration by parts) as well as negative r (by symmetry). 
Therefore, we have 

Y(a + ir, u) — 0(e -e ^) 

as r — > oo, uniformly in u. 

Second, let us consider the Dirichlet series associated with the set S , i.e., D(s ) = J2 h&s h~ s . It 
can be written as a product of elementary functions 

m 1 

dm =e =n 7—=;■ («) 

hes 3=1 


Each factor 1/(1 — Pj s ) has a simple pole at 0 and its singular expansion there is given by 
1 /( 1 ^Pj S ) ~ l/(slogPj) as s ->■ 0 . 
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Next we consider the Mellin transform of f(t,u) 7 which is given by Y(s,u) D(s). This function 
has a pole of order to + 1 at s = 0, so the Laurent series of Y(s, u)D(s) has the form 


fm(u) , fm-i{u) 


H-h 


fl(u) , fo(u) 


^ra+1 s ,,L S^“ S 

with f m {u) as indicated in the statement of our lemma. The other coefficients can be expressed in 
terms of certain improper integrals. Applying the Mellin inversion formula, we get 

(‘C+ioo 

+~ s ds 


i rc-\-i oo 

= — Y(s, u)D(s)t~ 

27r * Jc-ioo 


for any c > 0. Following Flajolet, Gourdon and Dumas |7 Theorem 4], we shift the line of integration 
to the left and pick up residues at the poles. This is possible because of the aforementioned growth 
properties of Y(s,u). The main contribution comes from the pole at s = 0, where the residue is 
indeed 

^Vgi/tr 1 

7 .' 


/l(zi)(logl/t) + /oO). 


, . _ . . 7 -T7y(l°g l/i) m 1 + - 

TO! (TO — 1)! 

There are further poles at all multiples of 2ni/logpj (1 < j < to), which are all simple poles (no 
two of them coincide) in view of the fact that the pj were assumed to be pairwise coprime, hence 
they only contribute 0(1). In fact, the 0(1) term can be replaced by a sum of to Fourier series 
with periods logp,,- (1 < j < to) that are given by 


\E 'j(t)= V"] Res Y(s 7 u)D(s)t 

* ^ s=2irik /log 

fcez\{o} ' 1 


E y ( 

fcez\{o} 


2nik 
log p. 


,u 


1 


n 


l°g Pj ^ 1 _ p - 2 -Gi° g P , 

r^j 


exp - 


/ 2irik log t \ 

V log/ ’ 


We remark that these Fourier series have exponentially decaying coefficients, since Y(s 7 u ) decays 
exponentially in imaginary direction, while Baker’s theorem on linear forms in logarithms (see 
Chapter 12 of 0) guarantees that 

m 1 

TT_L__ 

i-J- 1 1 -2-rrik/logpj i 
r= 1 | L ~ Pr 

r+2 

is bounded above by a power of k: indeed, there exist constants A^ 7 B a (depending on the bases 
Pi , P 2 j • • *, Pm) such that 

A = | k log p r — Clogpjl > A^k~ BA 

for all r ^ j and integers k and t not equal to 0. Thus, if || • || denotes the distance to the nearest 
integer, 


log p r 


log Pj 


and consequently 

|l _ ^-STri/c/logpj I 
| x rr | 

It follows that 


1 — expI — 


1 


> ~^k~ BA 

log Pj 


2mk log p r 
log Pj 


> 4 


klogPr 


II u —2irik/ logpj 
r=1 | 1 — Pr 


= O 


log Pj 


> k ~ BA 
log Pj 


which in turn means that 

Res Y(s 7 u)D(s)t~ s = o(\k\ {rn ~ 1)BA e~ 2 ^ kWlogp A . 

s=2irik/ logpj V / 


Thus each of the Fourier series is convergent and indeed represents a smooth function. This 
proves the asymptotic formula for f(t,u). The derivatives f(t,u) and J^/(f,w) have Mellin 
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transforms (1 — s)F(s — l,it)-D(s — 1) and (s — l)(s — 2)Y(s — 2,u)D(s — 2) respectively, so essentially 
the same arguments apply, now with the main terms coming from the poles at 1 and 2 respectively. 
It remains to prove the estimate for the third derivative. Note that it can be written as 


q3 

dt :i 


f(t, u)=Y h 


3 ^— ht 


Q(e~ ht ,u) 


h£S 


(1 


o ht . 


yd— lg— (d— l)ht ^3 


where Q is some polynomial. If we choose 77 (recall that our result will be valid for |Imt| < 77 ) 
small enough so that the denominator stays away from 0 (compare the analysis of Y(s,u) above), 
the last factor is uniformly bounded by a constant. The Mellin transform of 


J2 k 

h£S 


3 e~ ht 


is given by T(s)D(s — 3), to which we can apply the same arguments as for the harmonic sums 
encountered before. The dominant singularity is clearly a pole of order m at s = 3 in this case, so 
that the desired estimate follows immediately. 

The proofs of the second and the third part of Lemma are analogous. □ 

5. Estimating the Tails 

For our application of the saddle-point method, we need to estimate the tails (i.e., the parts 
where z is away from the positive real axis) of the generating functions given in (4.1 1 , (4.2 ) and (4.3). 
This is done in the following sequence of lemmas. First of all, let us introduce some notation. For 
r > 0 , we set 

S(r) = S fl [1,1/r] = {h £ S : hr < 1}. 

It is straightforward to prove that 

(log l/r) m 


|<S(r)| = 


■0(( logl/r)" 1 " 1 ) 


(5.1) 


m - nr=i 1 °spj 

as r —> 0 + . Note that later (starting with the next section), r will be determined by the saddle 
point equation. 

Lemma 2. Let u be in the interval [|,2], and let z = e ~ r+2niy with r > 0 and y £ \—\, |]- There 
exists an absolute constant C such that 

\FM\ 

F(\z\,u) 

\G(z,u)\ 


and 


G(\z\,u) 

I H b (z,u)\ 

H b { M 


<exp(-C INH' 

^ heS(r) 

<ex P f-c \\ h y\\ 

^ heS(r) 


^<exp(-C Y II^H 2 )- 


h£S(r) 

where || • || denotes the distance to the nearest integer. 

Proof. For positive real a and complex w, we have the two identities 


|1 + aw\ 

(1 + a |to|) : 


= 1 - 


2a(\w\ — Rein) 
(1 + a |w |) 2 


and 


1 + aw + aw 2 


(1 


= 1 - 


2a(|'u;| — Rew)(l + 2\w\ + a |w | 2 + 2Re w) 


( 1 - 


T 


Assuming that a £ [|,2] and |u;| < 2, we get 


|1 + aw\ 


(1 




< 1 - ^(h| - R ew) < exp(-^(|tw| - Re w)j 


(5.2) 
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and 


1 + aw + aw 2 


(1 + a |to| + a |w| ) 


9 < 1 — t —(\w 

2 '2 - 169 


-Rew) < expf—— (Iwl -Rew 
V 169 v 1 


(5.3) 


Now let d be even, set a = u and w = z h , so that (5.2 1 , together with the triangle inequality, yields 
1 + uz h + u 2 z 2h + • • • + u d - l z {d ~ l)h 

< |l + uz h \ + u 2 \z\ 2h |l + uz h \ H- 1 - u d ~ 2 I z \ < ' d ~ 2)h |l + uz h \ 


< (1 + u \z\ h + u 2 \z\ 2h + • • • + u d ~ l \z\ (d ~ 1)h ) exp(-^( \z\ h - Re(i 

Taking the product over all h G S gives 

\F(z,u)\ < F(\z\,u )exp(-l £ ( \z\ h - Be(z h ) )) 

u hes 

= F {\A , u) ex p(-^ Y e ~ hr ^ _ c °s(2t r/iy))) 


lies 


< F(|z| ,w)exp(-^- (1 - cos(2tt/i?/))) 

hGS(r) 

< F(|z|,w)exp(--^- Y INf), 

heS(r) 

which proves the first statement of the lemma with C = 4/(25e). For odd d , we can argue in a 
similar fashion, but we also apply (5.31 (with a = 1 and w = uz h ) and use the triangle inequality 
in the following way: 

\l +uz h +u 2 w 2 + ... +u d-l z (d-l)h\ 

< |1 + uz h + u 2 z 2h I + u 3 \z\ 3h |1 + uz h \ + ■ • • + u d ~ 2 U d ~ 2)h |1 + UZ h \. 


For the generating function G(z,u), the reasoning is fully analogous, but we also have to use (5.3) 
with a = u and w = z h . A similar situation occurs for Hb(z, u). □ 


Next we estimate the sum that occurs in the previous lemma. When m > 2, relatively simple 
estimates suffice for our purposes, while we need an additional auxiliary result in the case that 
to = 2. The following lemma provides the necessary estimates. 

Lemma 3. Let r > 0 and y € [— 5 , 5 ]. and set 

£ = E(r,y)= Y IIMl 2 > 

h£S(r) 

where again || • || denotes the distance to the nearest integer. For sufficiently small r, we have the 
following estimates for S. 

(a) If \y\ < r/2, then £ > Ai(y/r) 2 (log(l/r)) m ~ 1 for some positive constant A\ (that only 
depends on m and the set of bases {p±,p 2 , ■ ■ ■ ,Pm})- 

(b) If \y\ > r/2, then £ > A 2 (log(l/r)) m-1 for some positive constant A 2 (that also only 
depends on m and the set of bases {pi,P 2 , ■ ■ ■,Pm})- 

Now let to = 2. For any constant K > 0 and any 5 > 0, there exists a constant B > 0 depending 
on pi,P 2 , K and S such that the following holds for sufficiently small r. 

(c) We have £ > iFlog(l/r), except when y lies in a certain set E(K,r) of Lebesgue measure 
at most Br 3 ~ s . 


Proof. For better readability, the proof is split into several claims. 
A. Statement (a) is correct. 
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Proof o/[3} Let \y\ < r/2, which implies \hy\ < | for all h £ S( 7 '). Then we have 

s = II h yf = H h 2 y 2 > p 2 (y/ r f (\ s ( r )\-\s(r/p)\) 


h£S(r) 


hGS(r) 


hdS(r) 

h£S(r/p) 


for any p > 0. If we take p sufficiently small and apply the asymptotic formula in (5.1), we obtain 
estimate (a). □ 


B. A 2 can be chosen in such a way that statement (b) holds for \y\ < r 2 / 3 . 


Proof of\B[ Let us assume that r/2 < \y\ < r 2 / 3 . Then we have log 11 /y| > |log(l/r), and 
essentially the same idea as above works again. We obtain 

s = INII 2 > h2yl - H h2y2 > p 2 ( ,\ s ( 2 \y\)\~\ s (.\y\/p)\)’ 

h&S(r) heS(2\y\) fieS(2|y|) 

h ^s(\y\/p) 


and formula (5.1 1 can be applied again to obtain (b). 


□ 


We are left with the case that \y\ > r 2 / 3 , so we will assume this from now on. By Dirichlet’s 
approximation theorem, there exists a rational number a/q (with coprime a and q) such that 
q < r -2 / 3 and 



r 2/3 

< - . 

q 


C. There exists a positive constant Ci that only depends on m and the set of bases {pi,p- 2 , ■ ■ ■ ,Pm} 
such that for small enough r and any coprime integers a , q with q < r~ 2 / 3 , there are at least 

ci (log <?)(log l/r ) m_1 

many elements hi £ S(r x / 3 ) with ||a/ii/g|| > 1/q. 


Proof o/[^ For q = 1, the statement is trivial, so we assume that g / 1. Let us now distinguish 
whether q is in the set S or not. 

If q £ 6 >, then write q = p^pf 2 • ■ We have 


A = max(ai, a 2 , ■ ■ • ,a m ) > 


log q 

l0g(piJ?2 • ■ -Pm) ' 


Suppose that Oj = A. Consider the elements hi = p^p^ 2 ■ ■ -p^T £ S with 0 < /3j < o.i = A. For 
any of these hi, the number ahi/q is not an integer and thus ||a/ii/q , || > 1/q. Let us now find a lower 
bound for the number of such elements hi. Using (applied to the set Si = {s € S : p t \ s}), 
we find that for some positive constants ci and ci, there exist at least 


c 1 A|5 i (r 1 / 3 )| > dGoggXlogl/rp - 1 
elements hi £ S with hi < r -1 / 3 . 

If q S, then we clearly have || /j-io,/q , || > 1/q for all hi £ S, so the same statement as in the 
first case holds again. □ 


D. There exists a positive constant c that only depends on m and the set of bases {pi,p 2 ,... ,p m } 
such that for sufficiently small r and r 2 / 3 < \y\ < there are at least 

c(log l/r ) m_1 

many elements h £ S(r) with \\hy\\ > l/(3pi). 

Proof o/[SJ Let us divide the interval [1 /< 7 ,1/2] into subintervals 

Jo = [l/( 2 pi), 1 / 2 ], h = [ 1 /( 2 pi), l/( 2 pi)], ... 
whose ends have a ratio of pi (except possibly for the last one). There are at most 

log(g/ 2 )/log(pi) < c 2 log q 


such intervals. 
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By [C]and the pigeonhole principle, we can choose one of these intervals (Ij, say) such that for 
at least C 1 /C 2 (logl /r) m ~ l distinct numbers h\ £ S with hi < r -1 / 3 , the number ||/iia/g|| lies in 
this interval Ij, i.e., we have l/( 2 p{ +1 ) < ||/iia/<?|| < l/( 2 p}). 

Now we have 

hip^a 

q 

which means that we have at least C 1 /C 2 (log l/r ) m_1 elements h = hip{ £ S with ||o/i/g|| > l/(2pi) 
and 

h = hip{ < h\q < r -1 / 3 ?- -2 / 3 = 

r 

All of these numbers h are therefore in the set S(r). For sufficiently small r, it follows that 



hia 

= p\ 

q 


> 


2 pi’ 


\\hy\\ > 


ha 

q 


which proves the claim. 


r 2 / 3 h 1 r 2 / 3 hiq 1 , /, 1 

- >-— >-r 1/3 > -, 

q 2pi q 2pi 3pi 


□ 


E. A 2 can be chosen in such a way that statement (b) holds for \y\ > r 2 ^ 3 . 

Proof of\E\ The result follows from [D] since 

for A 2 = c/(9 pi) if r is sufficiently small. □ 


So (b) is now proven in both cases, and it remains to prove statement (c) of the lemma, so 
assume that m = 2. Choose some e £ (0, <5), set 

L= L(l-e)log Pl 1/rJ 

and define, for a positive integer M, the set 

D(M) = {u £ [0,1] : ||pfu|| < pf 2 for 0 < t < L with at most M exceptions}. 

The constant M will be chosen appropriately at the end of the proof. 

We get the following result, which almost proves (c). 

F. Set R = |_elog p2 1/rJ. If y is not contained in the set 

E = IJ {y £ [-|, |] : p 2 y mod 1 £ D(M)}, 

k<R 

then 

E > epf 2 M \og p2 l/r. 

Proof of [ 7 ] By our assumptions, there is no k < R such that p\y mod 1 £ D(M). Therefore, for a 
fixed k the inequality > pf 2 holds for more than M choices of l < L. Moreover, we have 

P\P 2 < r -1+e • r~ e = r~ 1 for all such k and £. 

It follows that 

2 = IIM 2 > Ji^yll 2 > {R+l)Mpf 2 > epf 2 M log P2 l/r, 

/iG<S(r) £<L k<R 

which is what we wanted to show. □ 


It remains to show that the set E is small. This is done in the following two claims. 

G. The Lebesgue measure of the set D(M ) is at most 0(^L M p^ J ^ L ). 
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Proof of~G\ First, note that ||pfa|| > pf 2 unless the (£ + l)-th and the (£ + 2)-th digit after the 
decimaP point in the pi-adic expansion of v are either both 0 or both p± — 1. For an upper bound, 
we relax this condition to both digits being equal. 

Therefore, for an element of D(M), at least L — M + 1 of the first L + 2 digits have to be equal 
to the previous digit. Allowing exactly j < M exceptions, there are ( L j^) number of ways to 
choose the “exceptional” digits. Moreover, each digit that has to be equal to the previous one 
reduces the Lebesgue measure by a factor of p\. 

Putting everything together, we end up finding that the Lebesgue measure of D{M) is at most 


which proves the claim. 



-(£+i)+j 
Pi 


0(1 


Mp M- L ) s 


□ 


We need one more claim, which concerns the size of the exceptional set E. 

H. The set E has Lebesgue measure C ) (r 1_e (log l/r) M+1 ). 

Proof o/|H| Since y £ [— |] (an interval of length 1) and p 2 is an integer, the Lebesgue measure 

A is preserved under taking the pre-image of v p 2 v mod 1. Therefore, we have 

\({y : p k 2 y mod 1 € D(M)}) = A(D(M)) 

and obtain 

A (E) < A (D(M)) = 0(RL M p™- L ) = 0(r 1_£ (log l/r) M+1 ) . 

k<R 

Note that the implied constant only depends on pi, p 2 , M and e. □ 


If we choose M = \Ke~ 1 p\ logp 2 l i then statement (c) follows from the claims above (in particular, 
[F and Hi with exceptional set E = E(K 7 r). Note that A (E) = 0(r 1_e (log l/r) M+1 ) = 0(r 1_<5 ). 
This completes the proof. □ 


6. Application of the Saddle-Point Method 

We are now ready to apply the saddle-point method (see Chapter VIII of p] for an excellent 
introduction), which gives us asymptotic formulas for the coefficients of the generating functions 
F(z,u ), G(z,u) and Hb(z,u). In the following, we use the notations ft(t,u ), ftt(t,u), ... for the 
derivatives of / with respect to the first coordinate. 

Lemma 4. Let u £ [|,2], and define r > 0 implicitly by the saddle-point equation 

n = - ft(r,u). 

The coefficients of F(z,u) satisfy the asymptotic formula 

[z n ]F{z, u) = 1 e ™~+^( r,M ) (l + 0((log n)"!" 1 - 1 )/ 5 )), 

\/27 xf tt {r,u) 

uniformly in u. Likewise, if we define r > 0 by 

n = - g t (r,u ), 

then the coefficients of G(z, u) satisfy the asymptotic formula 

[z n ]G{z,u) = 1 e rar+9(r ’ M) (l + 0((logn)- (m - 1 ) /5 )), 

V 2n 9tt(r, u) 

uniformly in u, and if we define r > 0 by 

n = -ha t t(r,u), 

2 We should rather correctly say “pi-point” instead of “decimal point” since p± is the base of our numeral system, 

but this may lead to even more confusion. 
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then the coefficients of H b fz,u) satisfy the asymptotic formula 

[z n ]H b { Z , u) = 1 e nr+h, b (r,u) (j + 0 ((jpg n ) ~ ( m ~ 1 )/ 5 ) ) . 

v 27r hb,a{r, U) 

Let us first give a short outline on the proof, which we only present for F, since the other two 
cases are analogous. We start by using Cauchy’s integral formula to extract the coefficient of z n 
from F(z,u). After the subsequent change to polar coordinates (z = we choose r to 

satisfy the saddle point equation. Thus the Taylor expansion in the central region simplifies (the 
first order term vanishes). Lemma [l] shows that r is of order (logn) m-1 /n. In the central region 
(to be defined later), the error term is 0(log(l/r) _ ( m_1 ^ 5 ) by Lemma [l] and we can complete the 
tails to get a Gaussian integral. The remaining parts of the integral are estimated by means of 
Lemmas [2] and [3] If m > 2, parts (a) and (b) of Lemma [3] already give sufficiently strong bounds. 
In the case that m = 2, we have to divide the tails further into a small “exceptional part”, where 
we apply (b), and the rest, where the stronger bound from (c) holds. 

So much for the overview; let us start with the actual proof now. 


Proof. By Cauchy’s integral formula, we have 

where C is a circle around 0 with radius less than 1. Let r > 0 and perform the change of variables 
z = e - * = e~^ r+ ’ lT \ so that the integral becomes 


i r 

[z n ]F(z, it) = — / exp(nr + f(r + it, u ) + inr) dr. 

27r J-k 


( 6 . 1 ) 


Now we choose r = r(n, u) > 0 to be the unique positive solution of the saddle-point equation 

n = -f t (r,u). (6.2) 

Let c be a constant such that (m — l)/3 < c < (in — l)/2, we choose specifically c = 2 (m — l)/5. 
Consider first the integral 


/o ~h 


rr( log 1/r) 

— r(log l/r)~ c 


exp (nr + f(r + it, u) + inr) dr. 


For |t| < r(logl/r) c , using Taylor expansion and Lemma [l] we have 

f(r + ir,u) = f(r,u) + if t (r, u)t - f tt (r,u) T - + 0[ | t | 3 sup | f tu (r + iy,u )|) 

Z ' l„l<T ' 


|y|<7 

= f{r,u) + if t {r,u)T-f tt (r,u) y + £>((log l/r) m ^ 1_3c ) . 


Therefore, by the definition of r in ( |6.2[ >, we have 

nr+f(r,u) /-r(log l/r)“° / J 1 


In = 


2ir 


r r( log 1/r) c / 2\ 

/ exp (-ftt{r,u)— ) dr (l + £>((log l/r) m_1_3c )). 

J— r(log l/r) -c \ ^ / 


Furthermore, 


/•r(log 1/r) c / t 2 

/ exp [-ftt(r,u)—) dr 

/ —r(log l/r) -c V ^ 

= / expf -ftt{r,u)— ) dr -2 


r(log 1 fr)~ 


exp( -ftt(r,u) — ) dr 


27r 


-2 


fttiy'i'U) Jr(log l/r)~ c 


ex P ( -ftt(r,u)— ) dr , 
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and 


0< [ exp (-f tt (r,u) T -) dr < [ exp(-^/ M (r, u)r(log 1/r) c ) dr 

J r(log 1 /r)~ c \ Z J J r(log l/r) -c K Z J 


(log 1/r)- 

2 exp(-/ tt (r, w)r 2 (log l/r)" 2c /2) 


ftt{r, u)r{\ogl/r)- 


= O 


(’’(log 1/r) 


-(m-l-c) - 7 (log 1/r 


-l-2c\ 


for a constant 7 > 0. Since m — 1 — 2c = (to — l)/5 > 0, the 0-term goes to zero faster than any 
power of log 1/r. Hence we have 

p nr+f(r,u) p nr+f(r,u) 

Io= . , -(l + Q((logl/r) m - 1 ~ 3c )) = ; (l + Offlogrp-t™- 1 )/ 5 )). (6.3) 


a/ 2 t r/ tt (r, u) 


\/2nftt(r, u) 


It remains to show that the rest of the integral in (6.11 is small compared to Jo- To this end, 
note for comparison that 1 /\J2tt u) is of order r(log l/r) - ^™ -1 - 1 / 2 . Now consider 


h = 


> r(log 1/r) - 


exp(w + f(r + ir, u ) + inr ) dr. 


Then 


l/i I <e nr+n ‘ 


’>«) f 

J r> 


r(log l/r) -c 


exp(Re(/(r + ir, u) — f(r, u ))) dr 


„nr+/(r,u) 


I F(e 


— (r+ir) 


> u )| 


F(e~ r ,u) 


dr. 


J r(log 1/r) 

If m > 3, then we can use Lemma[2]and parts (a) and (b) of Lemma[3]to show that the integrand 
\F{e~^ r+lT \u)\ /F(e~ r ,u) on the right hand side is d(exp(— CAi/(2it) 2 ■ (log l/r) m_1_2c )) for 
|t| < irr and d(exp(—CH2(logl/r) m_1 )) otherwise, which immediately shows that 


IT I = 


o(e” r+/ ( r ’“)(r exp (_C'H 1 /( 2 7i-) 2 • (log l/r)" 1 " 1 " 20 ) + exp(-C'H 2 (log l/r) r 


For to = 2, we need to be more careful. Again, part (a) of Lemma [3] can be used for the interval 
where |r| < tt r, with the same bound as above. The rest of the integral is split again: we choose a 
constant K > 0 such that CK > 1 (C as in Lemma [2]), and S > 0 such that S < CA 2 (A 2 as in 
Lemma [3]) . 

If y = —t/(27t) is not in the exceptional set E(K,r) as defined in Lemma|3j then we have 


|F(e 


— (r+ir) 


w)| 


F(e~ r , u) 


= 0(exp(—CK log 1/r)) = 0(r GK ) . 


By part (c) of Lemma [3j the set of r-values for which this estimate does not hold has Lebesgue 
measure C ) (r 1 ”' 5 ), and we have the estimate 


I F(, 


,-{r+ir) 


X )| 


= C>(exp(— CA 2 log 1/r)) = 0(r CA2 ) 


F(e~ r , u) 

for all those r. Putting everything together shows that 

\h\=0( e nr+f(r ’ u) (r exp (-CA X (log 1/r) 1/5 ) + r CK + r CA2+1 ~ s )), 

which again means that I\ is negligible, since the exponents CK and CA 2 + 1 — S are both > 1. 
The same reasoning can of course be applied to 


r»— r*(log 1/r) 


h = 


exp(?rr + f(r + it , u) + inr) dr. 


This finishes the proof for the function F(z, u). The proof for G(z, u) and Hb(z, u) is analogous. □ 
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7. The Number of Representations 
It is straightforward now to prove our main results. 

Proof of Theorems We specialize by it = 1 in Lemma ^ which gives us 

P{n) = \z n ]F(z, 1) = / e nro+/(r o , 1 )(i + 0 ((io gn ) — (m—D/ 5 ) ), 

where r 0 is given by the saddle-point equation n = —ft(r o, 1) (as its unique positive solution). 
Making use of Lemma [I] we get 


n = 


M 1) 

(m - l)!r 0 


(log l/r 0 ) m_1 + 0((log l/r 0 ) m_2 ) , 


which readily gives us 


log 1/ro = logn — (■m — 1) log log n — log 


■Mi) 

(m — 1 )! 


O 


f log log n 
\ logn 


(7.1) 


for n —>• oo. Now it follows that 


nr 0 = 


(log n) 


and Lemma [T] also yields 

f(ro, 1) = %^(logl/r 0 ) 


fra ( 1 ) 

(m — 1)! 

n . fm—1 ( 1 ) 


m—1 


1 + 0 


(log n) r 


to ! 

Mi) 

to ! 

/m—1 (1) 

(to — 1)! 


1 - 


(to — 1)! 
m(m — 1) 
log n 


(log l/ro) r 


log log n — 


/ log logn 
\ logn 


©((logn)™" 2 

Mi) 


log 


O 


f log log n 

V (logn) 2 


(log n) r 


1 + 0 


f log log n 
\ log n 


logn ° (to — 1)! 

©((logn) 171-2 log logn) . 


Since / m (l)/m! = n and / m -i(l)/(m — 1)! = kto(^” 1 1 logpj — logd)/2, this readily proves 
Theorem [I] (note that the factor / tt (r 0 ,1) only contributes O(logn) to logP(n)). 

To get the more precise formula (Theorem [n]) in the case to = 2, we only need to expand a little 
further. □ 

In principle, it would be possible to obtain similar (as in Theorem 0, more precise asymptotic 
formulas (in terms of log n and log log n) for all to > 2, but the expressions become very lengthy. 

8. Sum of Digits, Hamming Weight, Occurrences of a Digit 

This section is devoted to the central limit theorems for the sum of digits (Theorem |III| ) , the 
Hamming weight (Theorem |IV| ) and the occurrence of a fixed digit (Theorem [v| . We will only 
present the proof for the sum of digits; the other two proofs being analogous. The weak convergence 
to a Gaussian distribution will follow from the following general theorem (see [8j Theorem IX.13] 
and the comment thereafter, which states that it is sufficient to consider real values of u ): 

Lemma 5 (cf. [§] Theorem IX.13]). Let Xi,X 2 ,... be a sequence of discrete random variables 
that only take on non-negative integer values. Assume that, for u in a fixed interval H around 1, 
the probability generating function P n {u) of X n satisfies an asymptotic formula of the form 

Pn(u ) = exp(i?„(u))(l + o(l)) 

uniformly with respect to u, where each R n (u) is analytic in fh Assume also that the conditions 

r?vr ( \ 

R ' n { 1) + K{1) ^ °° aUd {R' n {l) + i?"(l))3/2 ^ 0 
hold uniformly in u. Then the normalised random variables 

X n - R' n ( 1) 

” {R'n{ 1)+R"(1)) 1/2 

converge in distribution to a standard Gaussian distribution. 
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Proof of Theorem \IU\ We use Lemma [4] Let X n be the sum of digits of a random multi-base 
representation of n, and let 

be the associated probability generating function. In the following, we write r{u) instead of just r 
to emphasize the dependence on u (of course, r depends on n as well). Moreover, we set Tq = r(l) 
as in the previous section. In view of Lemma [4] Lemma [5] applies with 

R n {u) = n(r(u) - r 0 ) + f{r(u),u) - f(r 0 , 1) - ^ log f tt (r{u),u) + ^ log f t t(r 0 , 1). 

We only have to confirm the conditions on the asymptotic behaviour of the derivatives. It is easy 
to extend the argument of Lemma [l] to obtain 

s ‘^( logl /i) m + 0(( logl /t) m -i), 


d j d k 

3W /(, -“ ) = 


du k 


3=0, 


( 8 . 1 ) 


.(-1 ) j (j - + 0(t^(logl/t) OT - 2 ) , j?0, 

as t — > 0+, uniformly in u. The definition of r by the implicit equation n = —ft(r(u),u) allows 
us to express r'{u) and all higher derivatives in terms of derivatives of / by means of implicit 
differentiation: we have r'{u) = — ftu{r{u),v) / ftt{r(u),u), and so forth. Thus it is possible to 
express the derivatives of R n only in terms of f{r(u),u) and its partial derivatives, for which we 
have the aforementioned asymptotic formula (18.11). Putting everything together, one obtains 


— R (, \ - — (— 
du k n to! \ du k ' 


1 


:/»(»)) (l»6 r( „) 

for k e {1, 2, 3}, so (making use of 0>) 

f'mi 1 ) 


o 




KW ~ ^(logl /roY 


(log n) m = 


and likewise 

<( 1) 


f" (1) 

^(logl/ro) r 


ml 


f"( 1) 

^(logn) m = 


-(u), 

1 

l°g Pi 

(d-l)(d-5) 


d — 1 

2m! 


n 


(log n) r 


ml 


12m! 


n i—(log n Y 

fJi log Pj 


and R'"{u) = 0((log?z) m ) uniformly in u. Thus the conditions of Lemma [5] are satisfied, which 
proves asymptotic normality of the distribution. However, we still need to verify the asymptotic 
behaviour of the moments (which is not implied by weak convergence). To this end, we apply the 
saddle point method once again. 

The generating function of the total sum of digits is F u (z, 1) = J^F(z,u) | , and the mean is 

given by 

[ z n }F u (z,l) 

^ [z n ]F(z, 1)’ 

so we have to determine an asymptotic formula for the coefficients of F u (z, 1). Cauchy’s integral 
formula, 

and the change of variables z = e _t = e~^ r ° +lT ' 1 (where ro satisfies the saddle point equation as 
before) yields 

i r 

[z n ]F u (z, 1) = — J exp(nr 0 +/(r 0 + ir, 1) + ott) /„(r 0 + zt, l)dr. 

Thus, 

[z n ]F u (z, 1) — /«(r 0 ,1) [z n }F(z, 1) 


i r 

= — y exp(nr 0 + f(r 0 + ir , 1) + inr) (/ u (r 0 + ir, 1) - f u (r 0 ,1 ))dt. 
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As we have seen in the proof of Lemma[4j the tails (the parts of the integral where |t| > r(log l/r) _c ) 
are negligible in that they only contribute an error term that goes faster to 0 than any power of 
log 1/r. So we may focus on the central part, where we expand into a power series 


exp(m~o + f{r 0 + it, 1) + inr) (f u (r 0 + it, 1) - / u (r 0 ,1)) = e nr ° +f(r °’ 1) Ft(r 0 d)T 2 /2 


fttui^Oi 1) 2 • 3 . ^ftttiXOi 0 F fttttuid (pi) 4 . 

T-r — i -r H-r + 

2 6 24 


We continue in the same way as in the proof of Lemma [4] to evaluate the integral over the central 
region asymptotically by making use of the asymptotic formula (8.1). This eventually gives us 

[z n }F u (z, 1) 


[z n Mz, 1) 
Thus in particular 


= fu(r 0l l) + 


ftu(r 0 , l)/ ttt (r 0 ,1) - ftt(r 0 , l)fttu(r 0 , 1) 


ftt(ro, l) 5 


O 


((log l/ r o) 


— (m—1) 


Pn = fu(r 0 , 1) + 0(1) = ^(log l/ro) r 

m\ 


O((logl/r 0 ) r 


c (d — 1) 


(logn)™ + 0((log?7,) m_1 loglogn) . 


21ogd 

We repeat the process with F uu (z, u ) + F u (z , u ) in the place of F u {z , u) to obtain an asymptotic 
formula for the second moment, which in turn yields formula 


2 [z n ](F uu (z, 1) + F u (z, 1)) 2 r / i\ i r / o , m(n i / \m—i\ 

V n = - [^n] F (-, 1) - Vn = fuu(r 0 , 1) + fu(r 0 , 1) + 0((log l/r 0 ) ) 


C(1) + /^(1) 

ml 

n{d — l)(d + 1) 


12 log d 

for the variance. This completes our proof. 


(log l/r 0 r + 0((log l/ro)" 1 " 1 ) 

(log n) m + O ((log n) 1 log log n) . 


□ 


9. Conclusion 


We obtained an asymptotic formula for the number of representations of an integer n in a 
multi-base system with given bases Pi, P 2 , ■ ■ ■, Pm, which are equivalent to partitions into elements 
of the set 

S = WP?---P'% n :aieNU{0}}. 

Moreover, we proved central limit theorems for three very natural parameters: the sum of digits 
(corresponding to the length of a partition), the Hamming weight (corresponding to the number 
of distinct parts of a partition), and the number of occurrences of a given digit. There are many 
more parameters that could be studied; to give one further example, the probablilty that the digit 
associated with a given element s € S in a random multi-base representation of n is equal to b for 
some b £ {0,1,..., d — 1} is 1/d in the limit as n —* oo, as one would heuristically expect. It is not 
difficult to adapt our saddle point approach to this problem, the generating function being 

1 - z hd 


^bs 


n 

hes 

h^s 


1 - z h 


in this case. As it was already mentioned in Section [2] it would also be possible to extend our 
results to other digit sets. 
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